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O ■ Abstract 



D , We solve the Einstein equations in the RandaU-Sundrum framework 

'^ ' with a static, sphericaUy symmetric matter distribution on the physical 

brane and obtain an approximate expression for the gravitational field 
^ ■ outside the source to the second order in the gravitational coupling. 

This expression when confined on the physical brane coincides with the 
standard form of the Schwarzschild metric. Therefore, the Randall- 
Sundrum scenario is consistent with the Mercury precession and other 
tests of General Relativity. 
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1. Introduction. 



Motivated by the hierarchy problems of the fundamental interactions, Randall and 
Sundrum have proposed an interesting scenario of extra non-compact dimensions in which 
four-dimensional gravity emerges as a low energy effective theory [1]. This novel proposal 
consists of two basic components. The first one is based on the assumption that ordinary 
matter and its gauge interactions are confined within a four dimensional hypersurface, 
referred to subsequently as the physical brane, embedded in a five-dimensional space [2]. 
The second ingredient involves the realization that a bound state of a five-dimensional 
graviton exists and is localized near the physical brane. Consequently, gravitational inter- 
actions at large distances on the physical brane are dominated by the bound state, and thus 
gravity appears effectively as four-dimensional. Modifications of this scenario include su- 
persymmetric generalizations [3] and models in which gravity is mediated with metastable 
graviton bound states [4] . Different issues concerning these models have been discussed in 
[5]. The Randall- Sundrum model consists of a three-brane embedded in a five-dimensional 
space which is asymptotically anti-de Sitter. A generic form of the metric in this space is 

ds^ = e-^^'^y^g^^dx^dx" + dy^, (1) 

the four-dimensional metric ^^j^ being asymptotically fiat. The asymptotic metric with 
g^u = ViJ.!^ satisfies the vacuum Einstein equations 

with the cosmological constant A = — 6k^. Here and throughout the paper, we adopt 
the convention that the Greek indices take values 0-3 and the Latin indices 0-4. The 
striking feature of this scenario is the existence of a massless graviton bound state which 
corresponds to the metric fiuctuations h^i, = g^i, — r]^^^, constant in y but nevertheless 
normalizable because of the conformal factor e"^**'^'. In the absence of nonzero modes 
(i.e., for a y independent metric g^v) the five-dimensional Einstein equations with negative 
cosmological constant reduce to the four-dimensional Einstein equations in vacuum without 
a cosmological constant. Thus the Randall- Sundrum scenario can be extended, provided 
that the metric g^i, is Ricci fiat, i.e., corresponds to any vacuum solution of General 
Relativity. This can be seen as follows: The Einstein equations can be rewritten as 

R^,. - 4«^(7/.. = 0, i?4^ = 0, i?44 - 4/^2 _ (3) 

with g^i, = e~^'^'y'g^iy. If we express the five dimensional Ricci tensor in terms of the four 
dimenional metric ^^j^, we find 

R^^u = R^,u + 'iK^g^^u■ (4) 
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for a y independent metric g^^. Consequently, g^i, obeys R^i^ = 0. An example is the 
Schwarzschild solution [6] 



ds^ = e-2'^l^l 



(1 )dt'^ + — 2GM ^^^ + ^^^^^ + ^y^^ (5) 



the source of which is a mass line extending infinitely in the y-direction. Another example is 
the pp-wave solution which describes gravitational waves propagating on the brane [7] . In 
the context of a physical situation-a spherically symmetric mass distribution located on the 
physical brane-nonzero modes are expected to contribute to the dynamics. The subsequent 
^/-dependence of the metric introduces difficulties in solving the Einstein equations exactly. 
How these modes localize the matter induced curvature near the brane yet maintain on 
the brane the asymptotic agreement between the solutions of the five-dimensional Einstein 
equations and the four-dimensional ones has only been shown explicitly for linear gravity. 
In this paper, we shall extend such an analysis to the first order of nonlinear ity ( second 
order in the gravitational coupling ) for the static axially symmetric metric produced by 
spherically symmetric matter that is localized on the brane. 

Similar to the case of linearized gravity, the solution to the Einstein equations can be 
formulated as a solution free from any metric singularity far away from the source when 
it is expressed in terms of coordinates that are straight with respect to the horizon. We 
shall eventually transform the solution to a new coordinate system-natural coordinates of 
the brane-by performing a gauge transformation that is determined by the Israel matching 
condition [8] ( Discussion on these issues concerning linearized gravity appear in [9]). We 
shall work in the weak field approximation and in a region outside the material source 
where the Israel matching condition reduces to the Neumann boundary condition. * 

Our solution, in contrast to the black-string solution [6] , will explore both the far away 
region on the brane and the AdS^ horizon. The arbitrary parameters of the solution are 
determined by the Newtonian limit on the brane. 

The apparent gauge dependence of the solution is the consequence of the incomplete- 
ness of the Randall-Sundrum scenario. Indeed, as it was pointed out in [10], all material 
particles fall freely towards the AdS^ horizon. Additional non-gravitational mechanisms 
are required to stabilize the brane and confine all material particles within the brane. In 
such models, possibly derived from string theory, the geodesic equation of material par- 
ticles would include extra contributions. But it is still plausible that the Israel matching 



* The physical brane contributes to the energy-momentum tensor as well. In the absence 
of matter, the metric ds^ = e~'^'^\y\ri^,ydx^dx'^ + dy"^ implies that T^^^ = mrifj,^6{y),T4^ = 
T44 = 0, where m is a constant. Its natural extention in the presence of matter is 
T^u = 'm-9fj,ryS{y)jT4^ = T44^ = 0. The linearity of Einstein's equation with respect to 
the second order derivatives leads to Neumann boundary condition. The same condition 
also guarantees the conservation of T^i^. 
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condition selects a gauge in which the motion of material particles within the brane is 
given by the simple four-dimensional geodesic equation 



In order to describe the real world, the Randall- Sundrum scenario has to satisfy all 
the existing tests of General Relativity. Within the linearized gravity, it has been shown 
that the Newtonian limit as well as the spatial part of the Schwartzschild metric can be 
recovered on the physical brane for nr » 1 with corrections from the nonzero graviton 
modes to the 1/r terms of the order 0{1/K^r^) [1], [11], [10], [12]. Thus the scenario 
implies a deflection angle of light within the experimental bounds. For the precession 
angle of Mercury, the observation accuracy is up to the first order of nonlinearity of the 
metric component goQ. Indeed, for the weak field expansion of the static spherical metric 



2GM G^M^ . ,.o ._ 2GM 



ds^ = -(1 7 ^^ + ...)dt'^ + (1 + + ...)dr^ + r^dn'^ (7) 



the precession angle measured in radians per revolution is given by [13] 

e- nGM , 

5 =-^(6 + 7) (8) 

where L is the semilatus rectum of the Mercury orbit and 7 = for the General Relativity. 
It is therefore imperative that the large r behavior of goo on the brane in the Randall- 
Sundrum scenario is examined beyond linearized approximation. Our perturbative solution 
reproduces the General Relativity result 7 = 0. Thus the Randall- Sundrum scenario is 
consistent with all experimental tests of General Relativity. 

For the sake of simplicity, we shall assume that the physical brane is located at y = 0. 
In the next section we shall derive the Einstein equations for a static axially symmetric 
five dimensional metric, investigate the gauge symmetries-coordinate trasformations that 
respect the particular form of the metric-and discuss the appropriate boundary conditions. 
In section three we shall obtain an exact solution to the linearized equations, a solution 
which is free of curvature singularities and describes the gravitational field produced by a 
spherically symmetric mass distribution located on the brane. In section four, we will derive 
an approximate expression of the metric to second order in the gravitational coupling (first 
order in non-linearity) which is valid in a region far from the mass source. This expression 
when confined on the physical brane is in agreement with the standard four dimensional 
Schwarzschild metric and thus reproduces the correct precession angle of Mercury. Finally 
we will discuss the physical implications of our results. 
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2. Einstein Equations, Symmetries and Boundary Conditions. 



The most general axially symmetric and static metric in D = 4 + 1 has the foUowing 



form 



-2K,\y\. 



ds^ = e-^'^i^l (-e^rft^ + e'dr' + e^r^dQ^) + dy' 



(9) 

where dfl^ = d6^ + sin 9d(l)'^ is the sohd angle on 5"^ and a, b and c are functions of r and 
y. Substituting the metric (9) into (3) we obtain the following equations 

-a" - -a' + -a' {-a' + b' - 2c') - e-^^y^^ \d - hnd -Kb- 2kc + -d(d + 6 + 2c)j = (10) 



a" + 2c" - -b' + -c' + -a' {a' - b') - c'{b' 
r r 2 



+ e 



-2Ky+b 



b — bub — Kd — 2k,c + -b{d + b + 2c) 



4 n' — h' 1 

c" + -c' + + c'2 + -{a'- b')c' 

r r 2 



+ e 



-2K.y+b 



c — K{d + b) — 6kc -\ — c(d + b + 2c) 



^b-c 



12 



d + 6 + 2c - 2K(d + 6 + 2c) + -(d^ + 6^ + 2c^) = 

a + 2c' - -(6 - c) + -a'(d - 6) - c'(6 - c) = 

r 2 



(11) 



(12) 

(13) 



where the prime denotes a partial derivative with respect to r and the dot denotes a partial 
derivative with respect to y. These equations apply to the positive side of the brane, i.e., 
y > 0, the corresponding equations to the negative side of the brane, y < 0, is obtained by 
switching the sign of k. 

For y-independent a, b and c, equations (13) are satisfied automatically and equations 
(10), (11) and (12) become those of a static, spherically symmetric and Ricci flat metric 
albeit not asymptotically anti-de Sitter. 



We shall discuss the gauge symmetries of these equations, more speciflcally the coor- 
dinate transformations that respect the axially symmetric, static form of the metric. Let's 
perform a coordinate transformation (gauge transformation) generated by v^u^ functions 
of rand y such that r i-^ r + v{r, y) and y ^^ y + u{r, y) and demand that this particular 
coordinate transformation respects the form of the five dimensional metric (9). We find 
that the functions v, u obey the following relations 



2{, + {,2 + ^-2Ky-2Kv+b{r+u,y+v)^2 ^ q 

v'{l + v)+ e-2'^^-2'^"+''('^+^'3'+")-u(l + w') = 



(14) 
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while the components of the metric transform under these residual coordinate transforma- 
tions as follows 

a{r, y) \-^ a{r + u,r + v) — 2kv 

b{r, y) ^b{r + U,y + v)- 2kV + In [(1 + u'f + ^2>^y+2KV-b{r+u,y+v)^>2^ ^^^^ 

u 
cir, y) I— > c(r + u,y + v) — 2kv + 21n (1 H — ). 

r 

In what follows, we shall construct the solutions to the Einstein equations in two different 
coordinate systems, or equivalently, in two different gauges [11], [14]. The solution in 
one coordinate system-the one in which the coordinates are straight with respect to the 
horizon-is free from metric singularities far away from the source. We shall denote this 
solution by a{r, y), b{r, y) and c{r, y) without superscripts and we shall impose the following 
conditions 

lim a = lim b = lim c = (16) 

\y\,r^oo \y\,r^oo \y\,r^oo 

independent of the order of the limits. The other set of coordinates are the natural coordi- 
nates of the brane. In this coordinate system the solution is given by the right hand side of 
(15) with the parameters u and v chosen appropriately to satisfy the Israel matching con- 
dition. We will denote the solution in this gauge by a^{r,y), b^{r,y) and c^{r,y), where 
the superscript P means "physical". Outside the material source, they satisfy the Neu- 
mann boundary condition, d^(r, 0) = b^{r,0) = c^(r, 0) = 0. In the subsequent sections, 
we shall use the following notation 

a^{r,y) = a{r,y) + 6a{r,y) 

b^{r,y) = b{r,y) + 6b{r,y) (17) 

c^{r,y) = c{r,y) +6c{r,y) 

where 6a, 5b and 5c are not necessarily infinitesimal. 

The Newtonian limit is specified by the asymptotic behavior of qqq on the brane for 
large r, i.e. 

a^r,y =1 18 

r 

with G being the four-dimensional gravitational constant and M the total mass of the 
sourse. For a weak gravitational field, GM « r, we have the following perturbative 
expansions for the solution that is free of metric singularities 



a{r, y) = ao(r, y) + ai(r, y) H , 6(r, y) = bo{r, y) + 6i(r, y) + ■ 

c{r, y) = co(r, y) + ci(r, y) H 

and for the physical solution 

a^{r, y) = a^ir, y) + af (r, y) H , 6^(r, y) = 6^(r, y) + 6f (r, y) + 

c^{r, y) = Co{r, y) + cf (r, y) H 



(19) 



(20) 
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Similarly we can expand the parameters of transformation 

v{r, y) = vo{r, y) + vi{r, y) + ■ ■ 
u{r, y) = uo{r, y) + ui{r, y) + ■ ■ 



f211 



All quantities with the subscript "0" are of the order GM and those with the subscript 
"1" are of the order G^M^ . In the next section we shall derive an exact solution to the 
first order in the gravitational coupling and we will discuss its physical significance. 

3. Linearized Solution. 



For a weak gravitational field, we obtain the linearized Einstein equations, i.e. 



- < - -/o - e-^'^l^l 


do — 5Kdo — tiiiQ — 2kco 


= 










a^' + 2c^'-^6;, + lc^ + e-2- 


bo — bubo — kcio — 2kco 


= 0, 








Co - K(do + bo) - 6kco 


-^{bo- 
r^ 


-co) = 


= 


ao + 6o + 2co - 2K(do + 6o + 2co) = 






2 • 
do + 2co (6o -Co) 


= 0. 

















(22) 



Initially, we shall seek solutions that satisfy equation (16). The fourth of these equa- 
tions is particularly simple, and the only solution consistent with the condition (16) is 



ao + bo + 2co = 0. 
Similarly, the last of equations (22) can be integrated to yield 



(23) 



a'o + 2co {bo- Co) = (j){r) 



(24) 



where (p is an arbitrary function of r. The boundary condition (16) then demands that 



Oq + 2co {bo -Co) = 0. 



(25) 



Note that the linearized Schwarzschild solution [6] satisfies (24) but not (25) . Substituting 
(23) into the first of equations (22) and using the Green's function found in [15], we obtain 



ao{r,y) 



8GM 2,„ p^ . . .K2{pe^y) 
-e^'y / dppjo{pr)- 



3nK 



Ki{p) 



(26) 
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where p = p/k. K^j{z) is the modified Bessel function of the second kind and jo{x) is the 
spherical Bessel function. The choice of the coefficient of the equation (26) is dictated by 
the requirement that it reproduces the correct Newtonian limit. Combining (23) and (25), 
we obtain a first order differential equation for h 

b', + ho + ^ = 0, (27) 

the solution of which reads 

It follows then from (23) then that 

, , 8GM ^^r^ ^ 1 . ^ ji{pr),K2{pe^y) ,,,, 

Co r,y = e^'^y / dpp-[jo{pr) - '-^^^] "f . 29 

37r«; Jo 2 pr Ki{p) 

It is straightforward to verify that the solutions (26), (28) and (29) satisfy the two of the 
remaining equations (22) which we have not used. 

If either r or C = -e'^^ becomes large, i.e, Kr » 1 or kz » 1, the integrals (26), 
(28) and (29) are dominated by the region where p << 1. The modified Bessel function in 
the denominator, Ki(p) ^ 4, and the integrals can be carried out explicitly. We find that 

4GM 2r' + 3C^ 4GM , 

-oir,0 = - 3 ^^,^^,^, =-^^(2 + cos e) 

, , ,, AGM 1 AGM 



2GM r^ + 2C^ 2GM ^ 

cor,C=^^— -^ ;^T^ = ^ — 1 + cos 9) 

3 (7.2 _^ ^2-) 2 3s 

where we have introduced four-dimensional spherical coordinates s and 9 via r = s sin 9 
and C, = s cos 9. In terms of s, and 9 the physical brane corresponds to cos 9 = —. We are 
primarily interested in the region with cos^ << 1. The result for ao(r, y) was obtained in 
[15]. Here we would like to emphasize that the expressions on the right hand sides of (30) 
approximate the entire region of ks >> 1. So it does interpolate the large r region on the 
brane and the large z region near the horizon. Furthermore, as is shown in the appendix, 
systematic inclusion of the higher order terms of Ki (p) will lend us a large s expansion of 
the type 

ig^y^ (31) 

n=0 

where Pn{9,\nKs) is a n-th order polynomial in In/ts with ^-dependent coefficients. 
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It is straightforward to check that while oq satisfies the Neumann boundary condition, 
bo and cq fail to do so. In order to transform the solution to the physical one we need to 
determine the parameters vq and uq. Equations (14) read to first order as 

vo = 0, v'o + e-^^yiio = (32) 

while the transformations of the components of the metric become 

Sao = —2kvoj dbo = —2kvo + 2uq, Scq = —2kvo + 2 — . (33) 

We can easily verify that the Einstein equations (22) remain invariant under the transfor- 
mations (32). Thus, the physical solution will be given by 

Oq = ao + Sao, bQ = bo + Sbo, Cq = co + Sco- (34) 

We impose Neumann boundary condition on b^ on the brane 6o'|y=o = 0. Under a coor- 
dinate transformation generated by vo,uo, we find that 

•p, -4GM ^ , -4GM ^ ,, ^ 

b^\y=o = ^-^ + 2< = -^^ - 2< = (35) 

Thus, 

Mr) = -i^—. uo{r,y) = —-^^e^^y + x{r) (36) 

6Kr bn'^r'^ 

where x is an arbitrary function of r and is chosen x = — 2C|M -^^ order to recover the 
standard Schwarzschild metric on the physical brane. The physical solution to first order 
takes then the form 

,P . !^e-. /" *pM!l) i^p^ + ?£M ^ |£^e-. (37) 

6T[K Jo pr J^iiP) or SK^r-^ 

6nK Jo 2 pr J^i[P) or Sk'^t'^ 

It is straightforward to verify that h^ = a^ + bQ + 2cq = 0. Using equations (30), the 
metric on the physical brane for r >> 1/k becomes 

ds^ = -(1 + ■ ■ ■)dt^ + (1 + + ■ • ■)dr^ + r2rfO^ (38) 

thus reproducing the standard form of the Schwarzschild metric. The dots in equation 
(38) represent terms of order 0{ '^ ^ ) and higher. 
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4. Second Order Approximate Solution. 



In this section, we will go beyond the first order solution (linearized) and find an ap- 
proximate solution to second order in the gravitational coupling of the Einstein equations. 
Inserting the ansatz (20) to equations (10)-(13) and keeping terms to second order, we find 



2,1, 



1 



a" a'l + -aQ(— ttg + &o — Scg) — e '^^ 'di — Bkcii — nhi — 2kci + -ao(tto + ^o + 2co) 



+ 6o(ao — 5Kdo — k6o ^ 2kco) 



ai + 6i + 2ci - 2K(di + 6i + 2ci) + -(d^ + 6o + 2co) 

2 • 1 

d'l + 2ci (6i - ci) + 7^ao(do - 6o) - Co(6o - co) = 



(39) 
(40) 







If we introduce hi = ai+bi + 2ci it follows then from (40) and taking into account equation 
(30) that 

kC'M^C [11 , 5 C^ ,3 C^ 



hi 



+ 



+ 



(^2 + ^2)2 Lis 9(^2 + ^2) 2(r2 + C2)2 

Consequently, integration over y yields 

G^M^ r25 7 C^ 1 C^ 



/ii- 



+ 



(41) 



(42) 



(r2 + C2)L36 18(r2 + C2) 4(r2 + C2)2. 
Taking into account (23), equation (40) becomes 

a'^h'^ - e-^^yb'oido - 4Kdo) + e-^'^^^/ii (43) 



a'l + -a[ + e-^'^y 
r 



di — 4Kdi 



with the right hand side completely known. We proceed to introduce spherical coordinates 

K 

takes the following form 



s, 6 such that r = scos^, z = -e'^^ = ssin^. In terms of these new variables, equation (40) 



1 a , 3 a, 1 

Ls"^ OS OS s^sm 



2.\r.'^ QQ0 



^ / • 2/1 ^ ^ 

-^""^a^^ sds ' .2 



3 a 3 , a 

+ — tant*-- 
oO 



ai{s,9) = fi + f2 + f3 (44) 



where 



2^^2 



/i = flo^o 



16^27^2 (2r2 + bC)r 



IQG'^M'^ 2 + cos2 9 - 3cos4 ^ 



9 



/2 = -e-2""6[,(do-4/.G 



(r2 + C2)4 



9 



KaQ) 



K' 



16G^M^ Scos^ 



(r2 + C 



2^4 



(45) 



/3 = e-2'^^/.(di+6i + 2ci) 



16G'2M2 1 
9 s^ 

10 



11 5 2., 27 4^ 

\ cos^ e H cos^ 

32 16 32 
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Equation (44) is a second order inhomogeneous linear differential equation, whose solution 
is a superposition of a special solution and a solution of the corresponding homogeneous 
equation. In order to derive the contribution to the solution from the inhomogeneous 
terms fi, i = 1, 2, 3 we make the ansatz (pi = ^i+Bi cos^j+c\ cos^ e ^ i = 1, 2, 3 with Ai, Bi, d 
constants. Then by demanding that 



4|-(^'|-) + ^^^(^^^'^^)--|- + 4tan^^l</>.(«,^) = /., z = l,2,3 (46) 
s"^ OS OS s^sm 9 OU OU s Os s^ 06 \ 

we calculate the coefficients Ai,Bi, Ci and thus determine the part of the solution which 
corresponds to the inhomogeneous terms 



?3 = 



G^M^ 8 4 2^ 8 4^, ^ 40^27^2 . 
^^(9 + 9^°^ ^+9^°^ ^)' <^2 = -y^^cos 9 

■(— + —cos"' 9 - 7Cos^ 9). 



(47) 



s2 '36 36 4 

Thus, a special solution of the equation takes the form 

<^ = </>! +</>2 +</>3 = ^(-| - ^cos^^ - ^COS^^). (48) 

s^ 36 36 12 

The most general form of the solution to the homogeneous equation, which is consistent 
with the large t^s approximation can be put in the form 

^[Fo{9) + F^{9) Inns + ... + Fm{9) In"^ ks]. 

As is shown in the appendix, the ordinary differential equation satisfied by Fm{9) can be 
converted into a hypergeometric equation with the two linearly independent solutions 

wi{9) = cos^6'cot6' 

and 

yj^(e) = --- ^cos^ 9 + -cos^ ^cot 9ln IZL^. (49). 

2 4 8 1 + sm 9 

wi is singular at 6* = (AdS^ horizon) and it should not be included while W2 is finite 
at ^ = ^, its inclusion though would produce on the brane terms of the form ^In"^ ^r, 
in disagreement with the standard Schwarzschild metric. Thus we expect F^{9) = 0. 
Similarly we expect that F^_i(6'), ..., Fi{9) are equal to zero and we are left with ai = 
(j) + r]W2, where 77 is a constant to be determined. 

Similarly to the solution of the linearized equations the physical solution to second 
order will be given by 

Oi = ai + 5ai, bi = bi + 5bi, c^ =ci + 5ci. (50) 

11 
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The parameters of transformation that generate the residual coordinate transformations 
vijUi obey the constraints (14) to second order 

2v, = -e-^^yiil = -e'^yv'i 

(51) 
t;; + e-2"^wi = -e-2'^^(6o - 2/ti'o)wo - e-2"^woWo 

where we have used equation (32), while the components of the metric transform to second 
order as 

5ai = OqUq + clqVq — 2kvi 

Sbi = bgUo + bovo - 2kvi + 2u[ - Uq + c'^'^^Vq 
dci = CqUq + cqVq - 2kvi +2 5- 



(52) 



and 

4 2?/^ 
5/ii = 5ai + bbx + 25ci = -%KVx + 2u\ + -wi - w[f + e^^^^v'^ ^. (53) 

The constraints (51) yield the following expressions for vi and iti, 

u^ = e^^yibo - 2tiVo + u'M - e^"^^! 

where A(r) is an arbitrary fucntion of r which can be determined by enforcing the brane 
bending condition- Neumann boundary condition to second order: more specifically, 



2 ITG^M^ 
— A — -j 



{h, + 5h,)\y=, = 2{-\" - -V - ^i^—^) = (55) 



We find that A(r) = — ^^^^2 and consequently that 

vi{r, y) = -e^'^y :—r ^ (56) 

A comment regarding the way Neumann boundary condition is met appears necessary at 
this point. Unlike the linearized solutions, which we have the exact expressions at our 
disposal, here only the approximate expressions at large ks are available. So we can not 
satisfy the Neumann boundary condition perfectly. On the other hand, it is easy to show 
that for large nr on the brane, di, 61 and ci are dominated by terms of order ^^4 . It 
is this order that has to be cancelled by the gauge transformation in conformity to the 
Neumann boundary condition. In another word, d^, b^ and c^ should be of higher order 
than — — :r- at y = and should vanish exactly when higher powers oi 1/ks are restored in 
the expressions of a^, b^ and c^ . 



12 
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Using equations (50), (52) and (56) we derive the following expression for af 

P G^M^ i25 16 ^ 17 2^ ,11 4 ^ 1 4 ^, 2 n 

tti = 7i — 1 CSC 6' CSC u + i 1 — CSC 6' esc yjcos 

^ s2 L36 9 36 ^36 3 18 ^ /57^ 

+ (— - -csc^)cos^ e + r]W2{e) 

It is straightforward to verify that af satisfies Neumann boundary condition d{'(r, 0) = 
only if ?7 = 0. On the physical brane and for nr » 1, the time component of the metric 
to second order reads 

-^00 = e"o +"1 =1 + 0( ^) (58) 

which is consistent with the General Relativity prediction 7 = 0. 

For completeness we will determine the expression for hi and ci on the brane and for 

Kr » 1. Using the fact that hi = ai + hi + 2ci, the second of equations (40) takes the 

form 

3 • • 1 • 1 1 

b[ + -bi = h[ + -hi di + -(ao'do + bo'bo + 2co'co) (59) 

r r r 2 

Furthermore, if we introduce (3 — bi — hi, equation (59) takes the simple form 

/3' + -P = fi + /2 (60) 

r 

where /i = — ^/ii — ^di and /2 = |(ao'do + boOo + 2co'co) Using now the new set of 
coordinates 8,9 and integrating (60), we find that 

(3 = ^ — ( cos^^- -cos^^)cos2^ (61) 

Inserting this expression into bi = P + hi and integrating, we derive 

bi = - —^ + — ^^cos^ 9. 62 

3 s'^ 12 s^ 

We still need to calculate the contribution from the residual coordinate transformations 
generated by vi,ui. According to equation (51), 61 transforms as follows 

bi^bi + b'^uo + boVQ - 2kvi + 2u[ - u'^ + e^'^^v'^ (63) 

with vi and ui given by equations (54). If we integrate the second of these equations, we 
derive an expression for ui 



4 ^2^2 1 G^M^ cot^ 9 
9 ssin^ 9 " 8 ~s sin^ 



^1 = a— ^7:- o^^-— TT+V' (64) 
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where t/j = t/j^r) is an arbitrary function of r. Combining equation (61) and the expression 
for Sbi, we find bf. The Neumann bondary condition is satisfied for arbitrary ifj. Imposing 
though the requirement that bf reproduces the standard Schwarzschild metric for large r 
determines t/^ 

43 G^ 

il^{r) = - — (65) 

72 r 

Combining relations (62), (63), (54), and (65), we find 



bf 



CM^n 8 . ^ 8 ^ 5 2^ 16 3^ 

— I — smc/H — cscc/H — esc esc 

-3993 9 



S2 



+ (— + -csc6' + — csc^ 6')cos2 e - -csc^ 6'cos'* 9 
^12 3 12 ^ 9 



(66) 



It follows from the definition of hf as well as the expressions of af and 6f , that 



c 






14 7 8 3 4 1 

1 — sin 6* csc^ 6 -\ — csc"^ 6 + ( — - H — sin 6* esc 6* 

6 9 6 9 ^43 9 

- — csc^ 9)cos^ 6* + (- + CSC 6* - — csc^ 6')cos^ 9 
12 ^ ^3 36 ^ . 



(67) 



What is remarkable here is that without any further degrees of freedom, we have already 

cf = 0(^^^) ^ (68) 

for cos^ = — and r >> 1. Thus, we have recovered the standard Schwarzschild metric at 

K,r ' 

large distances on the physical brane to the first order of nonlinearity 



2GM , . , 2GM 4G^M^ 
+ ■ ■ ■)dt^ + (1 + + , 



ds'^ ~ -(1 ^— + ■ ■ ■)df + (1 + ^^— + — ^ — + ■ ■ ■)dr'^ + r^rfO^ (eg) 



where the dots indicate terms of order 0(^^—3 — ) and higher. 
5. Conclusions. 



In this section we will recapitulate what we have done in this paper. We considered the 
most general ansatz for a static, axially symmetric metric in D = 4 + 1 dimensions in the 
Randall- Sundrum framework and derived the Einstein equations. We subsequently solved 
the equations to first order in the gravitational coupling and found an exact, singularity 
free solution of the linearized equations. This solution describes the gravitational field of 
a spherically symmetric mass distribution confined on the physical brane. At distances 
far away from the material source and on the physical brane this solution coincides with 
the four dimensional Schwarzschild metric. Thus we confirmed that all tests of Linearized 
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General Relativity are satisfied. Furthermore, we proceeded to examine the solution of 
the Einstein equations to the second order in the gravitational coupling GM. We derived 
an approximate expression for the gravitational field which is valid to second order in 
the gravitational coupling in a region far away from the source. When this expression is 
confined on the physical brane it coincides again with the standard Scwarzschild form to 
second order and thus reproduces the correct formula for the precession angle of Mercury. 
These results reinenforce the belief that the Randall- Sundrum scenario reproduces General 
Relativity on the physical brane to all orders in the gravitational coupling. We hope to 
report soon our results towards this direction. 
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Appendix 

The homogeneous equation (44), 

s'^ OS OS s^ sm 9 o6 o6 s os s^ 00 \ 

is a Quasi-Laplace since the first two terms inside the bracket is the four-dimensional 
Laplacian acting on an axially symmetric function. Inspired by the large s expansion in 
(31), we look for the solution of the form 



^ m 

-Y^Fm^r^Us (71) 



' 1=0 

with n a positive integer. Substituting it into (70), we obtain a system of m ordinary 
differential equations for F^'s 

^^7^g^{sm^e—) + 3ts.ne— + n{n + l)\Fi = il + l)[{2n + l)Fi+i-il + 2)Fi+2] (72) 

If we further introduce a new parameter g = cos^ 9, the differential equations (72) are 
transformed into 

g(l-g)^ + (-l-^g)± + ^n{n + l)]Fi = ^{l + l)[{2n + l)Fi+^-{l + 2)Fi+2] (73) 

Since Fm+i = Fm+2 = 0, the equation for Fm = F is a hypergeometric equation 

d'^F dF 

g(l-g) +(^-(a + (3 + l)g)—-a(3F = (74) 

dg^ dg 
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with a = ^, (3 = — i(n + l) and 7 = — 1. The other equations are of the same hypergeo- 
metric type but with inhomogeneous terms. The set of equations can be solved successively 
for Fi in a descending order of the subscript /. For odd n, /3 is a negative integer and there 
is a polynomial solution. The companion solution, however, is nonanalytic at ^ = or 
g = 1 or both. The terms of the large s expansion (31) contains odd powers of s and the 
power of the logarithm In ks does appear in association with the polynomial solution. For 
even n, the solution can be written as 

(75) 



with F satisfying a different hypergeometric equation with a = |(?t. — 1), (3 = —^ — 1 
and 7 = — 1. This new hypergeometric equation possesses a polynomial solution. But the 
square root in (75) does not make it an analytic solution for F. The case with 71 = 2 is of 
special interest. In this case one of the solutions is 

wi (g) = / = cos^ ^cot 9. (76) 

Given the solution wi we can derive the second solution by means of the standard formula, 

■-^ve/ -xj ^-^^i " ^ where p{ri) = — (i-^) • Subsequently we find that the 
second solution to the homogeneous differential equation is given by 

1 3 o^ 3 o„ „, 1 — sin6' 
= tCos"' 6 + -cos"* 6'cot 6'ln 

2 4 8 1 + sin^ 
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